Tissue-specific Transcriptome Research Articles

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects.

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

Nutri-cereal tissue-specific transcriptome atlas during development: Functional integration of gene expression to identify mineral uptake pathways in little millet (Panicum sumatrense).

Little millet (Panicum sumatrense Roth ex Roem. & Schult.) is an essential minor millet of southeast Asia and Africa's temperate and subtropical regions. The plant is stress-tolerant, has a short life cycle, and has a mineral-rich nutritional profile associated with unique health benefits. We report the developmental gene expression atlas of little millet (genotype JK-8) from ten tissues representing different stages of its life cycle, starting from seed germination and vegetative growth to panicle maturation. The developmental transcriptome atlas led to the identification of 342 827 transcripts. The BUSCO analysis and comparison with the transcriptomes of related species confirm that this study presents high-quality, in-depth coverage of the little millet transcriptome. In addition, the eFP browser generated here has a user-friendly interface, allowing interactive visualizations of tissue-specific gene expression. Using these data, we identified transcripts, the orthologs of which in Arabidopsis and rice are involved in nutrient acquisition, transport, and response pathways. The comparative analysis of the expression levels of these transcripts holds great potential for enhancing the mineral content in crops, particularly zinc and iron, to address the issue of "hidden hunger" and to attain nutritional security, making it a valuable asset for translational research.

Spinal cord injury-induced neurogenic bowel: A role for host-microbiome interactions in bowel pain and dysfunction

Background and aimsSpinal cord injury (SCI) affects roughly 300,000 Americans with 17,000 new cases added annually. In addition to paralysis, 60% of people with SCI develop neurogenic bowel (NB), a syndrome characterized by slow colonic transit, constipation, and chronic abdominal pain. The knowledge gap surrounding NB mechanisms after SCI means that interventions are primarily symptom-focused and largely ineffective. The goal of the present studies was to identify mechanism(s) that initiate and maintain NB after SCI as a critical first step in the development of evidence-based, novel therapeutic treatment options. MethodsFollowing spinal contusion injury at T9, we observed alterations in bowel structure and function reflecting key clinical features of NB. We then leveraged tissue-specific whole transcriptome analyses (RNAseq) and fecal 16S rRNA amplicon sequencing in combination with histological, molecular, and functional (Ca2+ imaging) approaches to identify potential mechanism(s) underlying the generation of the NB phenotype. ResultsIn agreement with prior reports focused on SCI-induced changes in the skin, we observed a rapid and persistent increase in expression of calcitonin gene-related peptide (CGRP) expression in the colon. This is suggestive of a neurogenic inflammation-like process engaged by antidromic activity of below-level primary afferents following SCI. CGRP has been shown to disrupt colon homeostasis and negatively affect peristalsis and colon function. As predicted, contusion SCI resulted in increased colonic transit time, expansion of lymphatic nodules, colonic structural and genomic damage, and disruption of the inner, sterile intestinal mucus layer corresponding to increased CGRP expression in the colon. Gut microbiome colonization significantly shifted over 28 days leading to the increase in Anaeroplasma, a pathogenic, gram-negative microbe. Moreover, colon specific vagal afferents and enteric neurons were hyperresponsive after SCI to different agonists including fecal supernatants. ConclusionsOur data suggest that SCI results in overexpression of colonic CGRP which could alter colon structure and function. Neurogenic inflammatory-like processes and gut microbiome dysbiosis can also sensitize vagal afferents, providing a mechanism for visceral pain despite the loss of normal sensation post-SCI. These data may shed light on novel therapeutic interventions targeting this process to prevent NB development in patients.

A multi-cell model for the C4 photosynthetic pathway in developing wheat grains based upon tissue-specific transcriptome data

A non-Kranz C4 photosynthesis of the NAD-ME subtype, specifically in developing wheat grains (14 dpa, days post-anthesis) was originally demonstrated using transcriptome-based RNA-seq. Here we present a re-examination of evidence for C4 photosynthesis in the developing grains of wheat and, more broadly, the Pooideae and an investigation of the evolutionary processes and implications. The expression profiles for the genes associated with C4 photosynthesis (C4- and C3-specific) were evaluated using published transcriptome data for the outer pericarp, inner pericarp, and endosperm tissues of the developing wheat grains. The expression of the C4-specific genes across these three tissues revealed the involvement of all three tissues in an orderly fashion to accomplish the non-Kranz NAD-ME-dependent C4 photosynthesis. Based on their expression levels in RPKM (reads per kilobase per million mapped reads) values, a model involving multiple cell- and tissue-types is proposed for C4 photosynthesis involved in the refixation of the respired CO2 from the endosperm tissues in the developing wheat grains. This multi-cell C4 model, proposed to involve more than two cell types, requires further biochemical validation.

Multi-omics resources for the Australian southern stuttering frog (Mixophyes australis) reveal assorted antimicrobial peptides

The number of genome-level resources for non-model species continues to rapidly expand. However, frog species remain underrepresented, with up to 90% of frog genera having no genomic or transcriptomic data. Here, we assemble the first genomic and transcriptomic resources for the recently described southern stuttering frog (Mixophyes australis). The southern stuttering frog is ground-dwelling, inhabiting naturally vegetated riverbanks in south-eastern Australia. Using PacBio HiFi long-read sequencing and Hi-C scaffolding, we generated a high-quality genome assembly, with a scaffold N50 of 369.3 Mb and 95.1% of the genome contained in twelve scaffolds. Using this assembly, we identified the mitochondrial genome, and assembled six tissue-specific transcriptomes. We also bioinformatically characterised novel sequences of two families of antimicrobial peptides (AMPs) in the southern stuttering frog, the cathelicidins and β-defensins. While traditional peptidomic approaches to peptide discovery have typically identified one or two AMPs in a frog species from skin secretions, our bioinformatic approach discovered 12 cathelicidins and two β-defensins that were expressed in a range of tissues. We investigated the novelty of the peptides and found diverse predicted activities. Our bioinformatic approach highlights the benefits of multi-omics resources in peptide discovery and contributes valuable genomic resources in an under-represented taxon.

S1 basic leucine zipper transcription factors shape plant architecture by controlling C/N partitioning to apical and lateral organs

Plants tightly control growth of their lateral organs, which led to the concept of apical dominance. However, outgrowth of the dormant lateral primordia is sensitive to the plant's nutritional status, resulting in an immense plasticity in plant architecture. While the impact of hormonal regulation on apical dominance is well characterized, the prime importance of sugar signaling to unleash lateral organ formation has just recently emerged. Here, we aimed to identify transcriptional regulators, which control the trade-off between growth of apical versus lateral organs. Making use of locally inducible gain-of-function as well as single and higher-order loss-of-function approaches of the sugar-responsive S1-basic-leucine-zipper (S1-bZIP) transcription factors, we disclosed their largely redundant function in establishing apical growth dominance. Consistently, comprehensive phenotypical and analytical studies of S1-bZIP mutants show a clear shift of sugar and organic nitrogen (N) allocation from apical to lateral organs, coinciding with strong lateral organ outgrowth. Tissue-specific transcriptomics reveal specific clade III SWEET sugar transporters, crucial for long-distance sugar transport to apical sinks and the glutaminase GLUTAMINE AMIDO-TRANSFERASE 1_2.1, involved in N homeostasis, as direct S1-bZIP targets, linking the architectural and metabolic mutant phenotypes to downstream gene regulation. Based on these results, we propose that S1-bZIPs control carbohydrate (C) partitioning from source leaves to apical organs and tune systemic N supply to restrict lateral organ formation by C/N depletion. Knowledge of the underlying mechanisms controlling plant C/N partitioning is of pivotal importance for breeding strategies to generate plants with desired architectural and nutritional characteristics.

Yolk proteins of the schistosomiasis vector snail Biomphalaria glabrata revealed by multi-omics analysis

Vitellogenesis is the most important process in animal reproduction, in which yolk proteins play a vital role. Among multiple yolk protein precursors, vitellogenin (Vtg) is a well-known major yolk protein (MYP) in most oviparous animals. However, the nature of MYP in the freshwater gastropod snail Biomphalaria glabrata remains elusive. In the current study, we applied bioinformatics, tissue-specific transcriptomics, ovotestis-targeted proteomics, and phylogenetics to investigate the large lipid transfer protein (LLTP) superfamily and ferritin-like family in B. glabrata. Four members of LLTP superfamily (BgVtg1, BgVtg2, BgApo1, and BgApo2), one yolk ferritin (Bg yolk ferritin), and four soma ferritins (Bg ferritin 1, 2, 3, and 4) were identified in B. glabrata genome. The proteomic analysis demonstrated that, among the putative yolk proteins, BgVtg1 was the yolk protein appearing in the highest amount in the ovotestis, followed by Bg yolk ferritin. RNAseq profile showed that the leading synthesis sites of BgVtg1 and Bg yolk ferritin are in the ovotestis (presumably follicle cells) and digestive gland, respectively. Phylogenetic analysis indicated that BgVtg1 is well clustered with Vtgs of other vertebrates and invertebrates. We conclude that, vitellogenin (BgVtg1), not yolk ferritin (Bg yolk ferritin), is the major yolk protein precursor in the schistosomiasis vector snail B. glabrata.

Transcriptomic profiling and the first spatial expression analysis of candidate genes in the salivary gland of the East Asian medicinal leech, Hirudo nipponia

Hirudo nipponia, a blood-sucking leech native to East Asia, possesses a rich repertoire of active ingredients in its saliva, showcasing significant medical potential due to its anticoagulant, anti-inflammatory, and antibacterial effects against human diseases. Despite previous studies on the transcriptomic and proteomic characteristics of leech saliva, which have identified medicinal compounds, our knowledge of tissue-specific transcriptomes and their spatial expression patterns remains incomplete. In this study, we conducted an extensive transcriptomic profiling of the salivary gland tissue in H. nipponia based on de novo assemblies of tissue-specific transcriptomes from the salivary gland, teeth, and general head region. Through gene ontology (GO) analysis and hierarchical clustering, we discovered a novel set of anti-coagulant factors—i.e., Hni-Antistasin, Hni-Ghilanten, Hni-Bdellin, Hni-Hirudin—as well as a previously unrecognized immune-related gene, Hni-GLIPR1 and uncharacterized salivary gland specific transcripts. By employing in situ hybridization, we provided the first visualization of gene expression sites within the salivary gland of H. nipponia. Our findings expand on our understanding of transcripts specifically expressed in the salivary gland of blood-sucking leeches, offering valuable resources for the exploration of previously unidentified substances with medicinal applications.

Genetic mechanisms and biological processes underlying host response to ophidiomycosis (snake fungal disease) inferred from tissue-specific transcriptome analyses.

Emerging infectious diseases in wildlife species caused by pathogenic fungi are of growing concern, yet crucial knowledge gaps remain for diseases with potentially large impacts. For example, there is detailed knowledge about host pathology and mechanisms underlying response for chytridiomycosis in amphibians and white-nose syndrome in bats, but such information is lacking for other more recently described fungal infections. One such disease is ophidiomycosis, caused by the fungus Ophidiomyces ophidiicola, which has been identified in many species of snakes, yet the biological mechanisms and molecular changes occurring during infection are unknown. To gain this information, we performed a controlled experimental infection in captive Prairie rattlesnakes (Crotalus viridis) with O. ophidiicola at two different temperatures: 20 and 26°C. We then compared liver, kidney, and skin transcriptomes to assess tissue-specific genetic responses to O. ophidiicola infection. Given previous histopathological studies and the fact that snakes are ectotherms, we expected highest fungal activity on skin and a significant impact of temperature on host response. Although we found fungal activity to be localized on skin, most of the differential gene expression occurred in internal tissues. Infected snakes at the lower temperature had the highest host mortality whereas two-thirds of the infected snakes at the higher temperature survived. Our results suggest that ophidiomycosis is likely a systemic disease with long-term effects on host response. Our analysis also identified candidate protein coding genes that are potentially involved in host response, providing genetic tools for studies of host response to ophidiomycosis in natural populations.

Tissue-Specific Transcriptomes in the Secondary Cell Wall Provide an Understanding of Stem Growth Enhancement in Solidago canadensis during Invasion.

Invasive plants generally present a significant enhancement in aboveground vegetative growth, which is mainly caused by variation in secondary cell wall (SCW) deposition and vascular tissue development. However, the coordination of the transcriptional regulators of SCW biosynthesis is complex, and a comprehensive regulation map has not yet been clarified at a transcriptional level to explain the invasive mechanism of S. canadensis. Here, RNA sequencing was performed in the phloem and xylem of two typical native (US01) and invasive (CN25) S. canadensis populations with different stem morphologies. A total of 296.14 million high-quality clean reads were generated; 438,605 transcripts and 156,968 unigenes were assembled; and 66,648 and 19,510 differential expression genes (DEGs) were identified in the phloem and xylem, respectively. Bioinformatics analysis indicated that the SCW transcriptional network was dramatically altered during the successful invasion of S.canadensis. Based on a comprehensive analysis of SCW deposition gene expression profiles, we revealed that the invasive population is dedicated to synthesizing cellulose and reducing lignification, leading to an SCW with high cellulose content and low lignin content. A hypothesis thus has been proposed to explain the enhanced stem growth of S. canadensis through the modification of the SCW composition.

Shifts in isoform usage underlie transcriptional differences in regulatory T cells in type 1 diabetes

Genome-wide association studies have identified numerous loci with allelic associations to Type 1 Diabetes (T1D) risk. Most disease-associated variants are enriched in regulatory sequences active in lymphoid cell types, suggesting that lymphocyte gene expression is altered in T1D. Here we assay gene expression between T1D cases and healthy controls in two autoimmunity-relevant lymphocyte cell types, memory CD4+/CD25+ regulatory T cells (Treg) and memory CD4+/CD25- T cells, using a splicing event-based approach to characterize tissue-specific transcriptomes. Limited differences in isoform usage between T1D cases and controls are observed in memory CD4+/CD25- T-cells. In Tregs, 402 genes demonstrate differences in isoform usage between cases and controls, particularly RNA recognition and splicing factor genes. Many of these genes are regulated by the variable inclusion of exons that can trigger nonsense mediated decay. Our results suggest that dysregulation of gene expression, through shifts in alternative splicing in Tregs, contributes to T1D pathophysiology.

Tissue-specific transcriptomes reveal potential mechanisms ofmicrobiome heterogeneity inan ancient fish.

The lake sturgeon (Acipenser fulvescens) is an ancient, octoploid fish faced with conservation challenges across its range in North America, but a lack of genomic resources has hindered molecular research in the species. To support such research, we created a transcriptomic database from 13 tissues: brain, esophagus, gill, head kidney, heart, white muscle, liver, glandular stomach, muscular stomach, anterior intestine, pyloric cecum, spiral valve and rectum. The transcriptomes for each tissue were sequenced and assembled individually from a mean of 98.3 million (±38.9 million SD) reads each. In addition, an overall transcriptome was assembled and annotated with all data used for each tissue-specific transcriptome. All assembled transcriptomes and their annotations were made publicly available as a scientific resource. The non-gut transcriptomes provide important resources for many research avenues. However, we focused our analysis on messenger ribonucleic acid (mRNA) observations in the gut because the gut represents a compartmentalized organ system with compartmentalized functions, and seven of the sequenced tissues were from each of these portions. These gut-specific analyses were used to probe evidence of microbiome regulation by studying heterogeneity in microbial genes and genera identified from mRNA annotations. Gene set enrichment analyses were used to reveal the presence of photoperiod and circadian-related transcripts in the pyloric cecum, which may support periodicity in lake sturgeon digestion. Similar analyses were used to identify different types of innate immune regulation across the gut, while analyses of unique transcripts annotated to microbes revealed heterogeneous genera and genes among different gut tissues. The present results provide a scientific resource and information about the mechanisms of compartmentalized function across gut tissues in a phylogenetically ancient vertebrate. Database URL: https://figshare.com/projects/Lake_Sturgeon_Transcriptomes/133143.

Brain and antennal transcriptomes of host ants reveal potential links between behaviour and the functioning of socially parasitic colonies.

Insect social parasites are characterized by exploiting the hosts' social behaviour. Why exactly hosts direct their caring behaviour towards these parasites and their offspring remains largely unstudied. One hypothesis is that hosts do not perceive their social environment as altered and accept the parasitic colony as their own. We used the ant Leptothorax acervorum, host of the dulotic, obligate social parasite Harpagoxenus sublaevis, to shed light on molecular mechanisms underlying behavioural exploitation by contrasting tissue-specific transcriptomes in young host workers. Host pupae were experimentally (re-)introduced into fragments of their original, another conspecific, heterospecific or parasitic colony. Brain and antennal mRNA was extracted and sequenced from adult ants after they had lived in the experimental colony for at least 50 days after eclosion. The resulting transcriptomes of L. acervorum revealed that ants were indeed affected by their social environment. Host brain transcriptomes were altered by the presence of social parasites, suggesting that the parasitic environment influences brain activity, which may be linked to behavioural changes. Transcriptional activity in the antennae changed most with the presence of unrelated individuals, regardless of whether they were conspecifics or parasites. This suggests early priming of odour perception, which was further supported by sensory perception of odour as an enriched function of differentially expressed genes. Furthermore, gene expression in the antennae, but not in the brain corresponded to ant worker behaviour before sampling. Our study demonstrated that the exploitation of social behaviours by brood parasites correlates with transcriptomic alterations in the central and peripheral nervous systems.

Genomic and transcriptomic analyses illuminate the molecular basis of the unique lifestyle of a tubeworm, Lamellibrachia satsuma.

Vestimentiferan tubeworms are representative members of deep-sea chemosynthetic ecosystems. In this study, we developed a draft genome and gene models and performed genomic and transcriptomic analyses of Lamellibrachia satsuma, the only vestimentiferan reported from the euphotic zone. The quality of the genome assembly and gene models is comparable to or higher than those of previously reported vestimentiferan tubeworms. Tissue-specific transcriptome sequencing revealed that Toll-like receptor genes and lineage-specific expanded bacteriolytic enzyme genes are highly expressed in the obturacular and vestimental regions, respectively, suggesting the importance of these tissues in defense against pathogens. On the other hand, globin subunit genes are expressed almost exclusively in the trunk region, supporting the hypothesis that the trophosome is the site of haemoglobin biosynthesis. Vestimentiferan-specific expanded gene families included chitinases, ion channels, and C-type lectins, suggesting the importance of these functions for vestimentiferans. C-type lectins in the trunk region, in particular, may be involved in recognition of pathogens, or in interactions between tubeworms and symbiotic bacteria. Our genomic and transcriptomic analyses enhance understanding of molecular mechanisms underlying the unique lifestyle of vestimentiferan tubeworms, particularly their obligate mutualism with chemosynthetic bacteria.

Tissue-specific transcriptome analyses in Drosophila provide novel insights into the mode of action of the insecticide spinosad and the function of its target, nAChRα6.

The insecticides spinosad and imidacloprid are neurotoxins with distinct modes of action. Both target nicotinic acetylcholine receptors (nAChRs), albeit different subunits. Spinosad is an allosteric modulator, that upon binding initiates endocytosis of its target, nAChRα6. Imidacloprid binding triggers excessive neuronal ion influx. Despite these differences, low-dose effects converge downstream in the precipitation of oxidative stress and neurodegeneration. Using RNA-sequencing, we compared the transcriptional signatures of spinosad and imidacloprid, at low-dose exposures. Both insecticides cause up-regulation of glutathione S-transferase and cytochrome P450 genes in the brain and down-regulation in the fat body, whereas reduced expression of immune-related genes is observed in both tissues. Spinosad shows unique impacts on genes involved in lysosomal function, protein folding, and reproduction. Co-expression analyses revealed little to no correlation between genes affected by spinosad and nAChRα6 expressing neurons, but a positive correlation with glial cell markers. We also detected and experimentally confirmed nAChRα6 expression in fat body cells and male germline cells. This led us to uncover lysosomal dysfunction in the fat body following spinosad exposure, and a fitness cost in spinosad-resistant (nAChRα6 null) males - oxidative stress in testes, and reduced fertility. Spinosad and imidacloprid share transcriptional perturbations in immunity-, energy homeostasis-, and oxidative stress-related genes. Low doses of other neurotoxic insecticides should be investigated for similar impacts. While target-site spinosad resistance mutation has evolved in the field, this may have a fitness cost. Our findings demonstrate the power of tissue-specific transcriptomics approach and the use of single-cell transcriptome data. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus

BackgroundIncreasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content.ResultsHere, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC.ConclusionsThe present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

Tissue-Specific Transcriptome and Metabolome Analysis Reveals the Response Mechanism of Brassica napus to Waterlogging Stress

During the growth period of rapeseed, if there is continuous rainfall, it will easily lead to waterlogging stress, which will seriously affect the growth of rapeseed. Currently, the mechanisms of rapeseed resistance to waterlogging stress are largely unknown. In this study, the rapeseed (Brassica napus) inbred lines G230 and G218 were identified as waterlogging-tolerant rapeseed and waterlogging-sensitive rapeseed, respectively, through a potted waterlogging stress simulation and field waterlogging stress experiments. After six days of waterlogging stress at the seedling stage, the degree of leaf aging and root damage of the waterlogging-tolerant rapeseed G230 were lower than those of the waterlogging-sensitive rapeseed G218. A physiological analysis showed that waterlogging stress significantly increased the contents of malondialdehyde, soluble sugar, and hydrogen peroxide in rape leaves and roots. The transcriptomic and metabolomic analysis showed that the differential genes and the differential metabolites of waterlogging-tolerant rapeseed G230 were mainly enriched in the metabolic pathways, biosynthesis of secondary metabolites, flavonoid biosynthesis, and vitamin B6 metabolism. Compared to G218, the expression levels of some genes associated with flavonoid biosynthesis and vitamin B metabolism were higher in G230, such as CHI, DRF, LDOX, PDX1.1, and PDX2. Furthermore, some metabolites involved in flavonoid biosynthesis and vitamin B6 metabolism, such as naringenin and epiafzelechin, were significantly up-regulated in leaves of G230, while pyridoxine phosphate was only significantly down-regulated in roots and leaves of G218. Furthermore, foliar spraying of vitamin B6 can effectively improve the tolerance to waterlogging of G218 in the short term. These results indicate that flavonoid biosynthesis and vitamin B6 metabolism pathways play a key role in the waterlogging tolerance and hypoxia stress resistance of Brassica napus and provide new insights for improving the waterlogging tolerance and cultivating waterlogging-tolerant rapeseed varieties.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

Tissue-specific comparative transcriptome analysis reveals candidate genes involved in lignans biosynthesis in Saururus chinensis

ObjectiveTo investigate the functional genes involved in the lignans biosynthesis in Saururus chinensis. MethodsAn RNA-Seq based tissue-specific comparative transcriptome analysis of S. chinensis leaves and rhizomes was performed, and an up-to-date bioinformatical analysis was carried out to discover the candidate functional genes. mRNA expression of part of the candidate genes were confirmed by qRT-PCR. ResultsAccording to the De novo transcriptome assembly, 49 430 unigenes in S. chinensis were presented and 33 222 unigenes (67.21%) were functionally annotated. Successively, 13 031 (26.36%) differentially expressed genes were filtered, and then GO and KEGG classification and enrichment analysis combined Swissprot function prediction revealed twenty-eight candidate genes encoding twelve enzymes involved in lignans biosynthesis pathways. The qRT-PCR results were consistent with the RNA-Seq analysis for the candidate genes. ConclusionThese results enrich our knowledge on lignans biosynthesis-related functional genes and facilitate further synthetic biology research on specific lignans in S. chinensis.